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Abstract 

We consider a possibility to naturally explain tiny neutrino masses without the lepton number 
violation. We study a simple model with SU{2)l singlet charged scalars (sf , s^) as well as singlet 
right-handed neutrino {vr). Yukawa interactions for Dirac neutrinos, which are forbidden at the 
tree level by a softly-broken Z2 symmetry, are induced at the one-loop level via the soft-breaking 
term in the scalar potential. Consequently neutrinos obtain small Dirac masses after the elec- 
troweak symmetry breaking. It is found that constrains from neutrino oscillation measurements 
and lepton flavor violation search results (especially for /i — )• 67) can be satisfied. We study the 
decay pattern of the singlet charged scalars, which could be tested at the LHC and the ILC. We 
discuss possible extensions also, e.g. to introduce dark matter candidate. 
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I. INTRODUCTION 



The neutrino oscillation data provide evidence that neutrinos have tiny masses 
which can be understood as a clear signature for physics beyond the standard model (SM). 
The simplest way to obtain neutrino masses may be to introduce S't/(2)j;^-singlet right- 
handed neutrinos, u^j {i = 1-3), which have Yukawa interaction with the SM Higgs boson. 
Then Dirac masses for neutrinos are generated after the electroweak symmetry breaking. 
However, in this naive mechanism, the Yukawa coupling constants for neutrinos have to be 
unnaturally small (< 10~^^) in comparison with those for the other fermions. The most 



familiar idea to solve the problem would be the seesaw mechanism {g] by introducing Majo- 
rana mass terms for z/jj. Taking the Majorana masses much larger than vacuum expectation 
value of the SM Higgs boson, very light Majorana neutrinos are obtained without excessively 
small Yukawa coupling constants. 

Another interesting possibility to avoid tiny Yukawa coupling constants is the radiative 
generation of Majorana masses for without Dirac mass terms. The original model was 
proposed by Zee [7|, in which Majorana neutrino masses are obtained at the one- loop level 
by the dynamics of the extended Higgs sector^. There have been several variant models 
in this direction Some of them ay , |l3| include dark matter candidates by 



imposing an unbroken Z2 symmetry which forbids the Dirac masses for neutrinos. In those 
models (seesaw and radiative ones), neutrinos are regarded as Majorana particles whose 
mass terms cause lepton number violating phenomena such as the neutrinoless double beta 
decay. Lepton number violation (LNV), however, has not yet been confirmed by experiments. 
Thus it is valuable to investigate a possibility that the tiny neutrino masses are generated 
in theories where the lepton number is conserved and Yukawa coupling constants are not 
extremely small. 

Radiative generation of masses for Dirac neutrinos would be an interesting possibility. 
There were several studies in past for such a scenario in various frameworks such as the left- 
right symmetry jl^. Iisl]. supersymmetry (SUSY) jl6|, and extended models within the SM 



^ If leptons couple with only one of two doublet scalar fields in the Zee model in order to eliminate the 
flavor changing neutral current interaction, the model cannot satisfy neutdno oscillation data 

^ In ref. |1(|, the second singlet fermion is added to the model in ref. |9| in order to satisfy neutrino 
oscillation data. 
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TABLE I: Particle contents of the ILDNM. Here L^, ^r, and $ are the S'?7(2)L-doublet fields of 
left-handed leptons, the right-handed charged leptons, and the 5[/(2)i-doublet scalar field in the 
SM, respectively. Three column from the right show particles added to the SM. 



gauge group [l7|, [iSj] (See also ref. [l9|). The simplest model seems to be the one in ref. |l7l |. 
where Dirac neutrino masses are generated at the one- loop level by introducing two SU (2)^- 
singlet charged scalar fields (sf , s^) as well as i'}^. In this letter, we show the one- loop Dirac 
neutrino mod el ( ILDNM) is compatible with neutrino oscillation data although this was 
overlooked in 17|]. We can find parameter sets which satisfy also constraints from searches 
for lepton flavor violation. We discuss the collider phenomenology of charged scalars under 
these parameter sets. Their decay pattern into leptons can be a characteristic feature of 
the ILDNM, by which the model could be tested at the LHC and the International Linear 
Colhder (ILC). 

We also discuss some extensions of the model briefly; accommodating dark matter can- 
didates, case with Majorana masses for z/^j, and so on. 



II. THE MODEL 

Particle contents of the ILDNM are listed in Table [B Three 5't/(2)i-singlet neutral 
fermions v\i {i = 1-3) are introduced such that Dirac masses for neutrinos exist. A softly- 
broken Z2 symmetry is imposed in the model so that Dirac masses can be forbidden at the 
tree level, where i/jj are assigned to be Z2-odd. Dirac neutrino masses are generated at the 
one- loop level by utilizing S'f/(2)i-singlet charged scalars, sf and s^, where is taken as 
a Z2-odd particle which can couple with The Yukawa interactions, the Higgs potential, 
and the Dirac neutrino mass generation in this model are presented below in order. 



3 



The Yukawa interactions for leptons are given by 



>CYukawa = ye Le ^ in + fw Ljicr2 Le + hei {inY u}^s^ + h.c, (1) 

where $ is the SM Higgs doublet field. The superscript c denotes the charge conjugation and 
CTj {i — 1-3) are the Pauli matrices. We take the basis where the Yukawa coupling matrix 
for charged leptons has been diagonahzed as y^. Notice that the matrix / is antisymmetric, 
while the matrix h takes somehow an arbitrary form. Although hu and fw are basically 
complex numbers, fw can be taken to be real numbers by using rephasing of three Li (and 

to keep real) without loss of generality. Furthermore we can take the basis where u}^ 
are mass eigenstates (of real positive mass eigenvalues). Then neutrino oscillation data give 
relations between the elements of fw and ha as shown later. 

The Higgs potential is written as 

V = v*^*+A(*w+/x^is^r+/x^i4r+{)"3 4«2 +h.c.} 

+ Ai + A2 141^ + {A3 {s+s^f + h.c.} + A4 14ri4l' 

+ A5($^$)|4r + A6($^<^>)l4l', (2) 

where //^ > 0. Although /x| and A3 can be complex parameters, they become real by 
rephasing sf and s^. Thus there is no complex parameter in the Higgs potential. Notice that 
/X3 is the soft-breaking parameter for the Z2 symmetry we imposed. The quartic coupling 
constants should satisfy the following conditions in order to avoid that the potential is 



unbounded from below; 

A > 0, Ai > 0, A2 > 0, (3) 

a;i = 2Ai + A5y^> 0, 2^^ + 2A3 + A4 + Ag^^ > 0, (4) 

uj2 = 2X2 + \e^>0, 2v'^ + 2A3 + A4 + A5y5> 0, (5) 

a;i2 = 2Ai + (2A3 + A4)y^>0, 2^^ + A5 + Aey^ > 0. (6) 



Mass eigenstates of two charged scalar fields are given by a mixing angle 6± as 

. -''''A (A , tan2.. = ^M^, (7) 
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FIG. 1: The one-loop diagram for Dirac neutrino masses in the ILDNM. Red dotted arrow shows 
the flow of lepton number. 



where we used m^± = fif + X-^v^ 12 and m^± ~ jl with v = a/2(0°) = 246 GeV. 

Clearly, > and > are necessary for (s^) = (s^) = 0. Masses of and H2 
are expressed as 



1 / / \ 2 



= - i "^l^ + - \l ( - ) + 4/^3 \ , (8) 



''1 



"^H± = \ \ "^L± + + \l ( "^!± - "^!± ) + 4/^3 }> > (9) 



where iff is defined as the lighter one. It is required to satisfy m ± m ± — /ig > so that 
rn? ± > at [s^) = (s^) = 0. The LEP experiment constrains masses of charged scalar 
fields to be greater than 73-107 GeV at the 95% confidence level (see mass limits for 



m 



from doublet fields and charged sleptons in ref. 

In Fig. [H the one-loop diagram^ for the Dirac neutrino mass is shown. The Dirac neutrino 
mass matrix M,^^ for {Mi,^)^ uf^ z/)j is obtained as follows: 

sin 20^ w± 

iM,,)u = Cif)ee'm,h,„ C^--^ln^. (10) 

^1 

Needless to say, 6± = and 7r/2 are not acceptable to obtain nonzero neutrino masses. Since 
we are taking the basis where i and z/)j are mass eigenstates, the Dirac mass matrix can be 



^ The same diagram was used in ref. 14] for the left-right symmetric model, where Yukawa coupling 
constant hu in the ILDNM is replaced by /«'. Current neutrino oscillation data cannot be satisfied in 
the case. 
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expressed as 



(11) 



where [i = 1-3) are neutrino mass eigenvalues {rrii > 0). The matrix [/mns is the Maki- 
Nakagawa-Sakata (MNS) matrix 21j which is expressed in the standard parametrisation 
as 
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-S23 C23 J 
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C12 S12 0^ 

-S12 C12 
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(12) 



where Qj and Sjj stand for cos 9ij and sin^jj, respectively. One of the m^'s is zero in this 
model because Det(M,^^) oc Det(/"^) = 0. Current neutrino oscillation data [l-5| allow two 
choices; either mi = or = 0. Two squared mass differences Amjj = mf — m?- are taken 



as Aml^ 



7.5 X lO'^eV^ and jAm^J 



2.3 X 10-3 eV^ 



For mi = 0, we have the following relations which were not found in ref. 17l |: 

f^e (f^MNs)^l 
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ITT^e (^MNs)e3 
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e3) 



(t^MNs)el 
(^MNs)^l 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



(f^MNSjel 

When we fix the MNS matrix and neutrino masses, six elements of hn and two elements of 
fu> in the left-hand side of these equations are given by five variables (/lei, he2, he3, f^r, and 
C). Rephasing of the massless u}^ makes hei real. Two phase degrees of freedom remain in 
hgi. However, eqs. (fT9|) and (l20l) show that the CP violating phase 6 in the MNS matrix 
vanishes because fw are real. 
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We assume for simplicity the so-called tribimaximal mixing 22| {s 



23 



1/2, 
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1/3,5^3 = 0) which agrees with neutrino oscillation data well. A simple example of the 
parameter set (the set A) which satisfy eqs. (fT3l) - (l20l) is 



h 



1.1 X 10" 



/ 



1.1 X 10" 



150 GeV, 



1 1.8 X 10"^ 1 

-2.4 X 10-^ -9.4 X 10"^ 1.0 X lO^^ 
1.4 X 10"^ -4.0 X 10"^ 5.9 X 10"^ 

0.5 0.5^ 
-0.5 1 
-0.5 -1 



(21) 



J 

200 GeV, ^± = 0.1rad. 



(22) 



(23) 



Notice that some elements of ha (especially, h^-i) tend to be small because of ratios of charged 
lepton masses in eqs. (fT3l)-(IT8l) while all elements of fw are in the same order of magnitude. 
For 1713 = 0, we obtain 
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(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



{UMNS)^i3 

Notice that eq. f l28p is the same as eq. (fT7|) . The phase of h^^ is absorbed by the massless 
z/|. while two phase degrees of freedom remain in /i^i and h^2- Equation ( 13T|1 means 6 = 
similarly to the case of mi = 0. It is worth to mention that we obtain = 1/2 and 5^3 = 
independently of hn if fe^ = —fer and f^r = 0, respectively. Such conditions on fu/ might 
be given by some discrete symmetry. Equations (|24|) -( 13T|1 for the tribimaximal mixing are 
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/i = 8.7 X 10 



-3 



satisfied by the following example (the set B) with values in eq. (123|) : 

( -7.0 X 10-1 1 ^ 

1 1 1 
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/o 1 -A 
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1 
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/ 



(32) 



(33) 



III. PHENOMENOLOGY 



In this section, we consider the constraint from the lepton flavor violating (LFV) decays 
of charged leptons and the prospect for the LHC physics. 



A. Lepton flavor violation 

The most stringent constraint on this model from the LFV decays of charged leptons is 
given by the experimental bound on the branching ratio (BR) of /i — )■ 67, BR(/i — > 67) < 
1.2 X 10^11 The branching ratio in this model is calculated as 

BR(/i ^ e7) ^ —V I [A- + -4-] + \ihh%ef(^ + -4-] \ '(34) 

7Q87TGI I \ml^ ml^J ^771%^ m^±; j 

where c± = cos6± and s± = sin^-|-. We ignore fermion masses in the loop integration and 
the electron mass in the final state. The parameter set A (eqs. ( l2T|) -( l23l) ) results in BR(/i — )■ 
67) = 2.9 X 10^12. This means not only that the ILDNM can satisfy the current bound on 
BR(/i — )■ 67) but also that the BR can be in the expected sensitivity of experiments in the 
future. On the other hand, the set B (eqs. f l23l) . fl32l) and fl33|) ) satisfies the experimental 
bound with a much smaller value BR(yU — > 67) = 7.5 x 10^^^. This is because f^r and /les 
for = are proportional to small S13. Although BR(r — fi'j) ~ 10"^^ for the set B is 
much larger than BR(/i — )■ 67), it is also far from experimental sensitivity. 

In the ILDNM the coupling constants and hu can be (9(10-^). Then the box diagram 
contributions to yU — )■ eee, which are proportional to the eighth power of these coupling 
constants, can be smaller than the current experimental upper bound although it becomes 
crucial in models where some of coupling constants are 0{1) {2^ . 
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B. Prospects at the LHC 



The charged scalar boson is expected to be produced at the LHC if it is light. 
The production cross section via qq — '~f*,Z* — )■ HiHi with ^/s = 14TeV is 23 fb for 
mjf± = 150 GeV for example. The partial decay widths of iu {£ = e, fi, r) are given 

by 

TTl zh / \ 

£' ^ \ e' i / 

where fermion masses are neglected. 

If is made dominantly from sf , its decay is determined by /«'. By using eqs. ( IT9ll 
and ( 120|) for mi = and the tribimaximal mixing, we obtain 

BR(iff -> ez/) : BR{H^ -> /ii^) : BR{H^ ^ rz/) = 2 : 5 : 5. (36) 

The set A (eqs. (!2T])- (!23l) ) gives approximately the same result. Equations (!30ll and ( 1311) for 
rris = and the tribimaximal mixing give 

BR{H{ ez/) : BR{H{ ^ ^z/) : BR(i/f ^ rz/) = 2 : 1 : 1. (37) 

The set B (eqs. (!32|) . ( 133|) . and (!23|) ) gives the same result in a good approximation. These 
results are robust because the matrix structure of f^^, is restricted very well. Therefore, if 
Hf" ~ sf, this model predicts BR(i7f — )• rz/)/BR(iff — )■ /iz/) ~ 1. On the other hand, if 
is made dominantly from sf, its partial decay widths are controlled by hii. For mi = 0, 
eqs. f ll3p and fll6p show h^-i ~ h^im^/m^. Furthermore, we have hrj, ~ h^'^m^/nij. with 
eq. f lTS]) for 6*13 = 0. For = 0, eq. fl2^ also means /i^3 ~ h^^m^/mr. Thus, it seems 
reasonable to expect hri ~ h^im^/mr {i = 1-3). Then we have 



BR(i7i" Tv) m 



2 

~ 3^ ~ 10-2. (38) 



BR(i7i /iz/) m 

As the result, this model is likely to give BR{H^ — )■ ru) /BR{Hi — )■ /iz/) < 1 due to the 
discussion above. 

If H2 is also light and the production cross section is significant, the decays into iu 
{£ = e,fi,T) smear the relation discussed above to some extent. Otherwise we can test the 
model at the LHC as well as the ILC by measuring the above characteristic pattern of the 
decay branching ratios. 



The partial decay width for h'^ — > 77 of the SM Higgs boson which is caused at 
the one-loop level in the SM, is affected by additional one-loop diagrams with and 
The contributions to the SM prediction depend on A5 and Xq as well as Higgs masses. 
When coupling constant for h^H^H^ (z = 1 or 2) is positive (negative), the additional loop 
effect from gives a destructive (constructive) contribution to the SM prediction. These 
contributions can amount to C(10) % deviations [25!]. If such the effect is detected at the 
LHC when the light SM Higgs boson is discovered, it can be an important indirect signature 
of the charged singlet scalar bosons. 



IV. DISCUSSIONS 



Possible extensions of the ILDNM are discussed in this section. First, we try to introduce 
dark matter candidates which do not exist in the model. Next, we consider the case with 
lepton number violation. Even if lepton number is not conserved and Majorana mass terms 
for z/)j are allowed, the mechanism to suppress the Dirac mass term is fruitful. 



A. Introducing dark matter candidates 

A possibility to accommodate dark matter candidates would be to impose an unbroken 
Z2 symmetry (we call it Z2) to this model in addition to the softly-broken Z2 symmetry 
such that all particles in the loop are Z2-odd'^. Since the SM charged leptons in the loop in 
Fig. [1] cannot be Z2-odd, they must be replaced by newly introduced Zg-odd fermions which 
can be understood as the fourth generation leptons. The Zg-odd Dirac neutrino could be 
the lightest Zg-odd particle which is stable. However, it cannot be identified as the dark 
matter because the spin-independent scattering cross section on a nucleon is too large to 



satisfy current data of direct searches [26|| due to the diagram mediated by the Z boson. 
Therefore, such a minimal extended model is excluded. 

We may consider the other model by taking different scalar particle contents, where the 
dark matter candidate enters and Dirac neutrino masses are induced radiatively. Such a 



Instead of the Z2 symmetry, lepton number can be used when it is conserved. For example, a fermion 
(boson) with a lepton number 2 (1) could be stable. From this point of view, lepton number conservation 
seems fit well for introducing dark matter candidates. 
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model can be found in ref. 18| where the following exact Zl^-odd particles are introduced: 

and a real neutral singlet scalar as well 
18| . a real scalar dark matter (Re(02) or s°) 



0\T\ 



an 5't/(2)i-doublet scalar field ($2 = (02 5 v^2 
as a neutral singlet Dirac fermion (N)- In ref. 



and the so-called Dirac leptogenesis 



27j via N decay are considered. If the Dirac fermion 



is the lightest Zg-odd particle, the dark matter is different from its anti-particle in contrast 
with the dark matter of the Majorana particle. It could be compatible with the asymmetric 



dark matter scenario 



The Dirac leptogenesis would be also achieved by the decay of 



Another simple possibility is the i?-parity-conserving SUSY extension where the lightest 
supersymmetric particle becomes a candidate for dark matter. The detailed study will be 
presented elsewhere. 

B. Baryogenesis 



There seem to be two possible extensions of this model in order to realise baryogene- 
sis, although the detailed analysis on them is beyond the scope of this letter. One is the 



electroweak baryogenesis 



29| . The scalar sector should be extended in order to have CP- 



violating phases and also to achieve strong first order phase transition at the electroweak 



symmetry breaking. The other is an application of the Dirac leptogenesis [27|]. It is known 
that the leptogenesis is possible without LNV. The number to be converted to the baryon 
asymmetry is free from the z/)j. number because the sphaleron does not act on the gauge 
singlet fields z/)j. 



C. Lepton number violation 



The mechanism to induce the Dirac mass terms for neutrinos can be applied also to the 



lepton number violating case, in which z/]j have Majorana masses as Mj(z/]j)'^z/]:j. Then, the 
type-I seesaw mechanism is realized at the two-loop level via the one-loop induced Dirac 



masses. This model was studied in refs. 
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By this loop suppression mechanism. 



Mi are much lighter than those in the tree-level seesaw mechanism. Consequently, Mj can 



^ A one-particle- irreducible two- loop diagram also exists for the Majorana masses of vl, which seems to 
be overlooked in ref. uj 
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be at the TeV scale without excessive fine tuning on couphng constants. Such TeV scale 
Majorana neutrinos could be tested at the LHC as well as the ILC. 

In the two-loop seesaw model, we can remove the soft-breaking term of the symmetry. 
Then the Majorana masses for vi^ are generated at the three-loop level and the lightest 
becomes a dark matter candidate. This model coincides with the model proposed by Krauss, 
Nasri and Trodden [9l. 



V. CONCLUSIONS 

We have investigated a simple model (ILDNM) with the mechanism for radiative genera- 
tion of Dirac neutrino masses without introducing lepton number violation. In the ILDNM, 
the Yukawa interaction L ^v^ is absent at the tree level because of the softly-broken Z2 
symmetry, so that it is induced at the one-loop level by the soft-breaking in the mixing be- 
tween and s^. Tiny neutrino masses are generated from the TeV scale dynamics. We have 
found the model can be compatible with the current neutrino oscillation data as well as LFV 
search results (especially for /i — )■ 67). There is no CP- violating phase in the MNS matrix 
in this model. It is possible that BR(yU — )■ 67) becomes large enough to be discovered by ex- 
periments in near future. The ILDNM is likely to give BR(iJf ~^ rz/)/BR(iJ]" — )■ ^iv) < 1. 
Characteristic features of are expected to be tested at the LHC and the ILC. We also 
have discussed several possible extensions of the model, which implement dark matter can- 
didates, mechanism for baryogenesis, and the radiative type-I seesawlike scenario by using 
one-loop suppressed Dirac masses. 
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